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We construct a family of self-similar Markov martingales with given marginal distributions. This 
construction uses the self-similarity and Markov property of a reference process to produce a 
family of Markov processes that possess the same marginal distributions as the original process. 
The resulting processes are also self-similar with the same exponent as the original process. 
They can be chosen to be martingales under certain conditions. In this paper, we present two 
approaches to this construction, the transition-randomising approach and the time-change ap¬ 
proach. We then compute the infinitesimal generators and obtain some path properties of the 
resulting processes. We also give some examples, including continuous Gaussian martingales 
as a generalization of Brownian motion, martingales of the squared Bessel process, stable Levy 
processes as well as an example of an artificial process having the marginals of CV for some sym¬ 
metric random variable V. At the end, we see how we can mimic certain Brownian martingales 
which are non-Markovian. 
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1. Introduction 

Constructing martingales with given marginal distributions has been an active area of 
research over the last decade (e.g. [1, 2, 9, 10, 14, 15]). (Here and in the entire paper, 
marginal distributions (also marginals) refer to the 1-dimensional distributions.) A con¬ 
dition for the existence of such martingales is given by Kellerer [12] (see Hirsch and 
Roynette [11] for a new and improved proof). 

Three constructions of Markov martingales with pre-specified marginal distributions 
were given by Madan and Yor [14], namely the Skorokhod embedding method, the time- 
changed Brownian motion and the continuous martingale approach pioneered by Dupire 
[4]. Recently, Hirsch et al. [10] gave six different methods for constructing martingales 
whose marginal distributions match those of a given family of probability measures. They 
also tackle the tedious task of finding sufficient conditions to ensure that the chosen family 
is indeed increasing in the convex order, or as they coined it, a peacock. 

In this paper, we deal with a different, albeit related, scenario. We do not start with a 
family of probability distributions, rather we start with a given martingale (the existence 
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of which is assumed) and produce a iarge famiiy (as opposed to just a handfui) of new 
martingaies having the same marginai distributions as the originai process. We say that 
these martingaies “mimic” the originai process. 

This same task was undertaken in [9] for the Brownian motion. It gave rise to the 
papers [1, 2] and [15] (who coined the term Faked Brownian Motion). Aibin [1] and 
Oieszkiewicz [15] answered the question of the existence of a continuous martingale with 
Brownian marginals. However, their constructions yield non-Markov processes. Baker et 
al. [2] then generalised Albin’s construction and produced a sequence of (non-Markov) 
martingales with Brownian marginals. In this paper, we extend the construction of [9] to 
a much larger class of processes, namely self-similar Markov martingales. 

Before formulating a solution to this problem we give a brief account on the origin and 
relevance of the mimicking question to finance, and more specifically to option pricing; 
that is the pricing of a contract that gives the holder the right to buy (or sell) the 
instrument (a stock) at a future time T for a specihed price K. The theoretical valuation 
of an option is performed in such a way as not to allow arbitrage opportunities - arbitrage 
occurs when riskless trading results in profit. The first fundamental theorem (e.g., [18], 
page 231) states that the absence of arbitrage in a market with stock price St, 0<t<T, 
is essentially equivalent to the existence of an equivalent probability measure under which 
the stock price is a martingale. (Here without loss of generality, we let the riskless interest 
rate be zero.) The second fundamental theorem (e.g., [18], page 232) implies that the 
arbitrage-free price of an option is given by E[(S't — Ar)+], where the expectation is 
taken under the equivalent martingale measure. In the classical model of Black and 
Scholes the stock price St is given, under the martingale measure, by the exponential 
Brownian motion St = Sgexp^aBt — ta^ 12). The parameter cr is known as the volatility, 
and is assumed to be a constant. The resulting expectation produces the well-known 
Black-Scholes formula for option prices. However, empirical evidence shows that in order 
to match the Black-Scholes formula to market prices of options one needs to vary cr. 
As a consequence it is natural to ask whether there exist alternative models of stock 
prices that result in a prescribed option pricing formula, such as Black-Scholes. Finally, 
it is easy to see that the collective knowledge of {E[(S't — K)~^],K > 0} determines the 
distribution of St or the marginal distribution, for example [8]. Therefore, if one wants 
to keep the option prices given by the original formula but without the limitations of 
the original process (such as constant volatility) one has to look for martingales (to have 
the model arbitrage free) with given marginals (to keep the same option prices). This 
question received much attention in the last ten years, see the pioneering work of Madan 
and Yor [14]. 

Throughout this paper, we assume that all processes are cadlag and progressively 
measurable. We will use the notation = to mean equal in distribution for random variables 
or equal in finite-dimensional distributions for processes, and this will be clear in the 
context. For a given random measure M{dx), the measure M(cdx) for c > 0 is defined 
by / g{x)M{cdx) = Jg{x/c)M{dx). We will also write E[M(da;)] to mean the measure 
defined by f ( 7 (a;)E[M(da;)] =E[J' g{x)M{dx)] for any positive function g. 

We start with a martingale Z which is also a Markov process, that is, for any bounded 
measurable function g, 'E[g[Zt)\Fs\ ='¥,[g{Zt)\Zs\ for s<t, where (J^ 4 )t>o denotes the 
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natural filtration of {Zt)t>o- We aim to construct new processes that will have the same 
marginal distributions as Z while retaining the martingale and Markov properties. We 
assume further that Z is self-similar, that is, there exists a (strictly) positive function q{c) 
such that {Zct)t>o = ((z(c)^t)t>o for all c> 0. For Z to be non-trivial and stochastically 
continuous at t = 0, we must have (see, e.g., [5]) that q{c) = for some k > 0, i.e., 

(^ct)t>o = (c"^t)t>o VoO 

(which implies that Zq = 0). We then say that Z is self-similar with exponent k, or 
simply, K-self-similar. 

Denote the transition fmrction of Z by P{s,t,x,dy) := V{Zt € dy\Zs = x) for s < t. 
Then the self-similarity of Z translates to the following scaling property on P: 

P{cs, ct, c'^x, dy) = P{s, t, x, dy) Vc> 0, s < t and x. 

If L is the infinitesimal generator associated with P, this scaling property is equivalent 
to 

CIT L qq'K f, —K — Ps; (P^) 

where tTc is the operator defined by iVcfix) = ficx). From this, we see that 

Lg — S 7Tg —K P \7TgK. . (1.2) 

In the following, we present a mimicking scheme to the process Z by randomising 
the transition functions. We will see that this is equivalent to time-changing the process 
together with an appropriate scaling. We then obtain some properties of the resulting 
processes and give some examples. 

2. Mimicking scheme 

Let Z be a Markov process with transition function P, which is self-similar with exponent 
K > 0. Note that if Z is a martingale, then it has a cadlag version; if Z is not a martingale, 
there are conditions for a Markov process to be cadlag. In this section, we construct new 
Markov processes that possess the same marginal distributions as Z and show that the 
resulting processes are martingales under certain conditions. 

Lemma 2.1. Let 0 < s <t <u. For any a,b G [0, 1] and a measurable set B, 
j P(0,s,0,da;)P(af,t, {at/s)'^x,B) = P(0,t,0,P) 


and 


P{at, t, {at/s^x, dy)P{bu, u, {bu/t)'^y, B) = P{abu, u, {abu/s)'^x, B). 
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Proof. Here we prove only the second equality, the first one is proved similarly. Suppose 
that a, & £ (0,1]. By the scaling property, we have 


P{at,t, {at/s)'^x,dy)P{bu,u, {bu/t)'^y, B) 


= / P{{bu/t)at, {hu/t)t, {bu/tY{at/sYx, {hu/tY dy)P{bu,u, {bu/tYy,B) 


= / P{abu,bu,{abu/sYx,dy)P{bu,u,y,B) = P{abu,u,{abu/sYx,B). 


Notice that when a or 6 is 0, the two equalities in the lemma trivially hold true. □ 


Proposition 2.1. Let {GsYs<t be a family of probability distributions on the set (0,1], 
where Gs,s = the Dirac measure at 1. Suppose that for any bounded measurable func¬ 
tion h and s <t <u, we have 

J y’h(a6)G,,t(da)Gt.„(d6)= y h(r)G,,„(dr). (2.1) 


Then P defined as follows is a transition function, 

PY>,t,0,dy) = P(f),tY,dy), 

P{s,t,x,dy) = J P{rt,t,{t/sYr'^x,dy)Gs,t{dr), s<t. 


Proof. Clearly, for each (s,t,x), P{s,t,x,dy) is a probability measure and for each 
{s,t,B), P{s,t,x,B) is measurable in x. Note also that P{s,s,x,B) = Sx{B). Next, using 
Lemma 2.1, we obtain, for 0 < s < t < u, / P(0, t, 0, dy)P{t, u, y, B) = P(0, u, 0, B) and 

J P{s,t,x,dy)P{t,u,y,B) = J J P{abu,u, {u/sY{abYx, B)Gs,t{da)Gt,u{db) 

= j P{ru,u,{u/sYr'^x,B)Gs,u{dr) = P{s,u,x,B); 

in other words, P satisfies the Chapman-Kolmogorov equations. □ 

Proposition 2.2. If Gs,t depends on s and t only through t/s (i.e. Gs,t =Gt/s), then 
the scaling property of P carries over to P: 

P{cs, ct, c'^x, dy) = P{s, t, x, dy) Vc> 0, s < t and x. 

Proof. This follows immediately from the definition of P and the scaling of P. □ 
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Let, for s <t, Rs,t be a random variable having distribution Gs,t- Property (2.1) is 
equivalent to the property that if, for s <t <u, Rs^t and Rt^u are independent random 

variables, then Rs,tRt.u = Rs,u- Further, if we let Va^b = — lni?ga gi, and write iFo,& for the 
distributions of Va,b (with Ka,a =<5o), then Property (2.1) translates to the convolution 
identity 

^a,b * — R-a^c^ a ^ b ^ C. 

As we seek to retain the scaling property of the original process, we assume that Gs,t = 
Gt/s and immediately reduce Property (2.1) to 

Ra * Rb — Ra+b^ U, 6 ^ 0. 

The family {Ra)a>o defines a subordinator (process with positive, independent and 
stationary increments) and by Levy-Khintchine it has Laplace transforms of the form 

J e~^'"Ra{dv) = exTp{—a'tp{X)), A > 0. 

The function ip, known as the Laplace exponent, takes the form 

poo 

V'(A)=/?A+ / (1 - e“'^'')z/(da::), (2.2) 

Jo 

with drift /3 > 0 and Levy measure i/ satisfying i^({0}) = 0 and /g°°(l A x)i^(da;) < oo. 
Conversely, to each ip (i.e., to each pair (/?,!/) as above) corresponds a convolution semi¬ 
group {Ra)a>o and in turn a family (Gu)u>i which satisfies Property (2.1). For details 
see, for example, [3], Section 1.2. This ensures the existence of (G„)„>i and a process 
with transition function P. 

Theorem 2.1. Let Z be a k- self-similar Markov process. To each ip, Laplace exponent 
of a subordinator, corresponds a k- self-similar Markov process X, starting from 0 and 
having the marginals of Z. Furthermore, if Z is a martingale and iP{k) = k, then X is 
also a martingale. 

Writing in terms of Rf/g, s <t, the new transition function 

P{s, t, X, dy) = W,[P{Rt/st, t, {t/sYR‘f/,.x, dy)], s < t 

can be seen as a randomisation of P{s,t,x,dy). Furthermore, the condition on {Gt/s)s<t 
for X to be a martingale can be written as = (s/t)^. 

Proof. By the Kolmogorov extension theorem, there exists a Markov process X with 
transition function P{s,t,x,dy). As for the martingale property, we first observe that 

^(A)=-iln /e-^’'K,(dn) = --ln /r^Ge.(dr), 
a J a J 
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so that V'(^) = ^ translates to J HGu{dr) = 1 for m > 1. Then we have 
J yP{s,t,x,dy) = J {t/s)'^r'^xGt/sidr) =x 

using the definition of P and the martingale property oi Z. □ 


The process X can also be obtained using subordination (with a suitable scaling in the 
state space and an appropriate time-change). The idea of subordination was suggested 
by Bertoin^ in the context of Brownian marginals. However, subordination alone (albeit 
with a suitable scaling in the state space) is not sufficient. A further logarithmic change 
of time is needed. This naturally creates an issue at 0. To deal with this, we follow 
Oleszkiewicz [15]. 


Proposition 2.3. Let (Ct)t>o be a subordinator independent of Z. Let, for a G K, 




.^a+ln 


t > e““. 


Then the process (A(‘^^)j>e-a has the same marginal distributions as {Zt)t>e-’^ and there 

exists a process {Xt)t>o such that for any a € R, iXt)t>e-‘‘ = (^i'^^)t>e-“ • Tbe process 
(Xt) t>o is a K-self-similar Markov process with transition function 


Q{0,t,0,dy) = P{0,t,0,dy), 

Q{s,t,x,dy) = (t/s)”e“”^'"<*/'*>a;,d?/)], s<t. 

Moreover, X is a martingale provided that Z is a martingale and = (s/t)'^. 


Proof. Since Zt{ui) is measurable in uj and right-continuous in t, it is measurable as a 
function of {t,uj). Hence, for each t > e““, x[°‘\uj) is a random variable. 

For c> 0, let Z^ = e“'‘^‘=Z^eC= and C* = Cc-i-s - Cc so that (.?s)s>o = iZs)s>o and 
(Cs)s>o = (Cs)s>o- Then we have, for t > e~^ > e““ and with c = a — b, 


(xn 


= (t'‘e“'‘f'’+'"*Z - '' _Arv(b)' 


t>e 


e^fc+lnt h>e 


)t>e-^ = (xn 


t>e 


For e““ < ti < • • • < let be the law of (A(“\ ... ,Xj“''). Then the measures 

(/iti,...,t„)n,o<ti<'"<t„ are consistent and by the Kolmogorov extension theorem, there 
exists a process {Xt)t>o with finite-dimensional distributions {pti,...,tn)n,o<ti<---<tn- 


^(a). 


A similar argument shows that {X^‘f^,.. .,X^'l) = 


,A. 


from which we 


deduce that (At)i>o is K-self-similar. As such, X extends by continuity to t > 0 by letting 


^Private communication. 
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Xq = 0. The equality of marginal distributions of and Z follows from the scaling 

property of Z as = t'^Zi = Zt for any fixed t > e-“. 

Using successively Lemmas A.3, A.4, A.2 and A.l (see the Appendix), we see that 
XM is Markovian. By the scaling property of P and the stationarity of subordinator 
we obtain, for e““ < s <t, the transition function of X^°''> as 

(s, t, X, dy) = e d 2 /|X(“) = x) = E[P(e"^‘“(‘/“U, t, dj/)]. 

Notice that the transition function does not depend on a and it is the same 
as P defined earlier with Gt/s being the distribution of exp(—The rest of the 
assertions of the proposition then follows immediately from Theorem 2.1. (Alternatively, 
we can carry out the proof independently, without referring to P, following Oleszkiewicz 
[15].) □ 

The process constructed in Proposition 2.3 is identical (in law) to the process ob¬ 
tained in Theorem 2.1 with Gt/g being the distribution of exp(—or R^/g = 

exp(-Cin(i/^*))- 

Remark 2.1. For Ai to be a martingale, /3 in (2.2) must be at most 1, and is 1 if and 
only if t (zz = 0 and '0(A) = A for any A) and X = Z. 

Remark 2.2. In Proposition 2.3, we cannot replace Ca-i-int with a two-sided subordinator 
Ct = Ct lt>o — C-tlt<o for t G R, where (Ct^)t>o and {Cf)t>o are independent subordinators. 
This is because by doing that, we will not have independent increments. In particular, 
since Co = 0, then for f < 0, the increment Co ~ C* = ~Ct G Gt, where Qt denotes the 
filtration generated by (. 


3. Properties 

In this section, we obtain the infinitesimal generators of the process X and display some 
of their path properties. We will work within the martingale framework, that is, unless 
otherwise stated, we will assume that our initial process Z is a martingale and we will use 
a subordinator C with drift j3, Levy measure v and Lapace exponent satisfying 0 (k) = n. 

Theorem 3.1. Suppose that Z has infinitesimal generator Lt- Then the infinitesimal 
generator of the process X is given by 

Aq/(x) = Lof{x), 

At fix) = PLtfix) -k (1 - /3)|x/'(x) 

+ 7 / [ ifiy)-fix))Pite~'^,t,xe~'^'^,dy)n{du), t>0, 

^ J (0,oo) J —CO 

for f differentiable and in the domain of L. 
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Proof. First, from Lemma A.3, Zt := e~*'^Zet is time-homogeneous with generator 
Lf{x) =Lif{x) — Kxf'{x) and transition semigroup Ptf{x) = //(?/)P(e“*, 1,dy). 
Next, applying Lemma A.4, the generator of the process Zt-=Z(;^=e is 

Lf{x)=l3Lif{x)-KPxf'{x)+ [ /'(/(y)-/(x))P(e"“,l,xe"“'',dy)z/(du). 

J (0,oo) J 


Then, let Zt = e'^^* Zt = e'^i* “^e and using Lemma A.2, the generator of Z is 

Ltf{x) = /37rg-,.(t-a)Pi7rg,.(t-a)/(x) + (1 - /3)ra/'(x) 

+ / /(/(e"(‘-“V) - f{x))P{e-\ 1, e-'^(‘-“)xe-“", dy)i/(du), 

J (0,oo) J 

since 7rg-„(t-a)A7rg,.(t-a)/(x) = A/(x) for Af{x) =xf'{x). Finally, we time-change the 
process Zt with a + lnt to get x[°'\ Thus, by Lemma A.l, the generator of Xt°'^ is 

= PLtfix) -b (1 - /3)|x/'(x) 

+ 7 / f ifiy)- f{x))P{te~'^,t,xe~'^''^,dy)v{du) 

^ J{0,oo) J 

due to a change of variable, the scaling property of P and identity (1.2). The generator 
of Xt is established by noting that does not depend on a. 

Since Z is self-similar, tTc/ is in the domain of L for all c > 0 whenever / is in the 
domain of L. Therefore, / is in the domain of A, if / is also differentiable. □ 

Note that when /3 = 1 and i/ = 0, we recover the process Z and At = Lt- 
For a measurable function /, if there exists a measurable function g such that for each 
t, fg |y(^s)| ds < oo almost surely and f{Xt) — f{Xo) — Jg g{Xs) ds is a martingale, then 
/ is said to belong to the domain of the extended inhnitesimal generator of X and the 
extended infinitesimal generator Asf{Xs) = g{Xs). If f{x) = x^ belongs to the domain 
of the extended infinitesimal generator of X, then X has predictable quadratic variation 

{X,X)t= [ AJ{Xs)ds. (3.1) 

Jo 

See examples in Section 4 for the computation {X, X) in some specific cases. 

Proposition 3.1. Suppose that Z is continuous in probability. Then the process X is 
also continuous in probability, that is, for every t, 

VoO, limP(|Xt-Xs| >c)=0. 
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Proof. We have, for s,t>e-\ 

F{\Xt -Xs\>c)< - s"e-"«“+'-)^e«.+.n* I > 

However, the first term 


(3.2) 


p(^|(t"e-"«“+'"* - I > 

< P(^|(t« - s'^)e-<“+'-Z I > fj 

+ p(^|s"(e-<“+'"‘ -e-<“+-)^e«.+i„J > l), 

which converges to 0 as s —> t, since f is continuous in probability as a subordinator. 

To deal with the last term in (3.2) we first observe that a process that is continuous in 
probability does not jump at fixed points so that P(Ae^“+‘"* = 0) = 1. Further, since Z 
is also continuous in probability, lims^t P(|.^yt — ZyJ > e) = 0 as soon as lims_>.t ys = yt- 
Therefore, for s < t + 1, 




= E 


liinPj |ZeC„+,„j 1) l^t ) l|Ae^‘*+ln‘=( 


=0} 


= 0 , 


where = o'(Cs, s < a + In t). 


□ 


Proposition 3.2. If Z is continuous in probability with finite second moments and 
has no drift (/3 = 0), then X is a purely discontinuous martingale. 

Proof. Let Ut = e^“+‘"* so that Xt = t'^Uff'^Zjj^ for t > e““. First, we observe that with 
probability one, Z does not jump at Ut- if U jumps at t. Indeed, if F is a countable set 
of points in [0,oo), then, as Z is continuous in probability, P(3t G F s.t. XZt 0) = 0. 
Let A = {t > e““: AH* > 0} and Fj/ = U~{A) where Uf = Ut-, then 

P(3t > e-“ s.t. AUt > 0, AZ^u,-) 7^0)= E[P(3s € F^ s.t. AZs ^ Ol^oo)] = 0, 

where Q denotes the filtration of U. Taking a to infinity, we obtain the desired result. 
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To show that X is purely discontinuous, that is, = 0, we compute the sum 

of the square of jumps of X. In general, we have 


E 

E 

= E 

E iAXg)'^lAUs>0,AZ^u,_)^0 


■e““<s<t 


■e““<s<t 


+ E 


-E 




At/s >0, 


‘e 


E (AXs)^1ac/,=o 


As Z is continuous in probability, the first term is zero due to the observation at the 
start of the proof. We write l{a,t) = for the third term. 

For the second term, we have, on the set {Alls > 0, = 0}, 

(AA,)2 = - UrA)f 

= s^^{U-^^{Zu^ - Z^u._)f + - U:^f 

+ 2ur{ur - U7_^)Z^u._){Zu^ - Z(c;,_))). 

Let 9t = E[Z^]. Since Z is K-self-similar, 6t = . As Z is a martingale, M[Zs {Zt — Zg)] = 

0 and E[(Zt — Zg)^] = ^[Z^] — E[Zg]. Thus, we obtain 

E[(AA,)2|aoo n {AC7, > 0, AZ(c/^_) = 0}] 

= - 0(o._)) + 0^u._){Ur - U;^f) 

= 2s2«6»i(l-C/7'=t/;_). 


Since {AZ([/^_) = 0} has probability one on the set [AUg > 0}, it follows that 


E 


E {AXg)'^lAu, 


>0,AZ([j^_)=0 


e““<s<t 


= E 


= E 


s>0 


E 2s2-0i(i_c/7"C7;_)1ac/. 

6““<s<i 

^ 2e2«(’-fo0i(l-e-"^«'-)lAc.>o 


^ 0 <r<a-\-lnt 

Writing in terms of v and tp, and using (2.2) with ip^n) = k, 

^o+ln t 


E 


_ ^o-|-inr p 

E {AXgflAUMz^u^_,=o =0i 2e2''(’'-“)dr / (l-e-^")Kd2) 

^ . ^0 «/(0,cx)) 


6““<S<t 


= 6>i(t^'^-e-^“'=)(l-/3). 
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Adding those three terms and taking limit as a —>■ oo, we have 


E 


E 


^0<s<t 


= — P) + lim l{a,t). 

a—foo 


Since X is square integrable on any finite interval if Z has finite second moments, it has 
quadratic variation with expectation E[[A, X]*] = E[Xj^] = \ Therefore, 


E[(A^A=)J=E[[X,A]i]-E 


E 

0<s<t,AXs#0 


= lim l{a,t). 


Note that both E[(X‘^, and l{a,t) are non-negative. Thus, when /3 = 0, we must 
have hmQ_).oo lia,t) = 0, which gives E[(A'^, A'^)t] = 0 and hence {X'^,X'^)t =0. □ 


4. Examples 

Given any self-similar Markov martingale Z with transition function P and infinitesimal 
generator L, we can mimic Z as per Section 2. We construct a new Markov martingale 
X that has the same marginal distributions as Z and possesses the same self-similarity 
Z enjoys from each f. We assume that the subordinator C has drift /3 and Levy measure 
V with Laplace exponent ip satisfying iP{k) = k, or = (s/t)'^. 

Examples of subordinators include Poisson process, compound Poisson process with 
positive jumps, gamma process and stable subordinators. For example, we can take C 
as a Poisson process with rate k/( 1 — exp(— k)) to satisfy iP{k) = k. In the following, we 
provide some examples of mimicking with the infinitesimal generators and the predictable 
quadratic variations computed explicitly to have a better understanding of the processes. 

We finish this section with a discussion on modifying our construction to mimic some 
Brownian related martingales and its limitation. 

4.1. Gaussian continuous martingales 

For any fc > 0, define the process {Zt)t>o by Zt = /J di?s, where {Bt)t>o is a Brownian 
motion. Note that 

iZt)t>o = (-B(i/(2fc+i))t2fc+i)t>o. 

This is a Gaussian process with zero mean and covariance function Cov(Zt, = 

^ Markov process and also a martingale. Moreover, it is (fc -I- ^)-self- 
similar {k = k since for all c > 0, 

{Zct)t>o = (5(i/(2fe-i-i))c2'=+it2'=+i)t>o = (c^"'’^^^-B(i/(2fc+i))t2fc+i)t>o = 
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A key aspect of the construction in [9] is the following representation of the mimic Xt 
when k = 0: 

At = t>s, 

where i?s,t has distribution Gj/s, ^s,t is standard normal and, Rs,t, ^s,t and Xg are 
independent. This representation extends to the case of other Gaussian continuous mar¬ 
tingales. In fact, it also extends to the case of stable processes - see Proposition 4.4. 

Proposition 4.1. With K = k + ^, the mimic {Xt)t>o has the representation 

Xt = {t/sY{RltXs + s"(l - Rs,tTis,t), t > s > 0, 

where Rs,t = , ^s,t = Zi and, Rs,t, ^s,t and Xg are independent. 


Proof. Since Z and C have independent and stationary increments, for e “ < s < t, 

where C is an independent copy of C, and ^g^t is a random variable distributed as Zi. 
Note that this representation holds also for t = s. □ 

Knowing that Z has generator Ltf{x) = \t^^f"{x) for t > 0, we can compute the 
generator of X following Theorem 3.1 and obtain, for t > 0, 

Atfix) = -pt^^nx) + (1 - P)^^xf{x) 

+ 7 / r ifiy) - /(x))P(te-“,t,xe-('=+i/2)-,dy)j.(du). 

^ J (0,oo) J — oo 

Taking f{x) = x^, then Equation (3.1), along with routine calculations and Equation 
(2.2), gives the following result. 


Proposition 4.2. The predictable quadratic variation of X is 

{X,X)t = + + 1) + (2fc + 1 - i,{2k + 1)) ds. 

Since {Z, Z)t = fg s^^ ds, we can also write 

(X,X}t = —^f^(2k + l)(Z,Z}t + (2k + l-^(2k + l)) f ixfds. 

ZK + i Jo S 


Remark 4-1- When k = 0, Z is a Brownian motion and our results agree with [9]. 
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4.2. Martingale of squared Bessel process 

A process {St)t>o is a squared Bessel process of dimension S, for some (5 > 0, if it satisfies 
dS”; = 2^/StdBt + (5dt, where B denotes a Brownian motion. The squared Bessel process 
S started at 0 is a continuous Markov process satisfying the self-similarity with k = 1. 
Let Zt = St — St. Then {Zt)t>o is a 1-self-similar Markov process and satisfies the SDE 


dZt — 2\/Zt St dBt- 


Note that S is stochastically dominated by the square of the norm of an n-dimensional 
Brownian motion, where n> S, thus E[5't] < nt and £[(2, Z)t] < 2nt‘^. It follows that Z 
is a true martingale (see, e.g., [13], Theorem 7.35). 

The infinitesimal generator of Z is Ltf{x) = 2{x St)f"{x), t>0, thus, that of X is 


Aof{x) = Lof{x) = 2xf"{x), 

Atf{x) = 2P{x + Sf)f"{x) + ^(1 - IS)xf'{x) 



“,d?/)j^(du), t>0. 


Proposition 4.3. The predictable quadratic variation of X is 



t>0 


where ip is the Laplace exponent of f. 


Proof. From dZt = 2\/Zt -I- St dBt, we have for w > 0, 



'2 

te~ 



te~ 


Zg \/Zg -|- Ss dBg 4 


Zgds + 2St^{l-e-^^) 


and taking conditional expectation (jJ Zg\JZg SsdBg is a true martingale - see above). 


E[Z^\Zt,-.] = Zl-^ + 4t(l - e-“)Zte-. + 2St^{l - e-"“). 
Thus, using Equation (2.2) we obtain, with /(x) = x^. 


Atf{x) = 2Stip{2) 4x(V'(2) - V'(l)) + ^^^(2 - V'(2))- 


However '0(1) = 1. The result then follows from Equation (3.1). Note that if /3 = 1, 
0(A) = A for any A and we recover (X, X)t = 2St^ -TdJ^Xg ds. □ 
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4.3. Stable processes with 1 < ck < 2 


Suppose {Zt)t>o is an a-stable process with 1 < a < 2. Then Z is a Markov process and 

{Zc’^t)t>o = {cZt)t>o, Vc> 0, 

that is, Z is K-self-similar with k = It is a Levy process with Levy triplet (0, 1 ^ 2 , 7 ), 
where 

vz{dz) = (4l^>o + dz 

for some positive constants A and B. Assume that Z is a martingale, in which case the 
Levy triplet must satisfy 


7 + 



zvz{dz) = 0 . 


Proposition 4.4. The mimic {Xt)t>o has the representation 

Xt = {t/sT{Rl^X, + s"(l - Rs,tT^s,t), t > s, 

where Rs,t = , ^s,t = Zi and, Rs,t, ^s,t and Xg are independent. 

Proof. See Proposition 4.1. □ 

The stable process Z has infinitesimal generator 

Lf{x) = + [ {f{x + y)- fix) - yfix)lDiy))nz{dy), 

Jr\{o} 

where D = {x: |x|<l}. To distinguish the characteristics of ( from that of Z, we add 
the subscript C to the drift and Levy measure of the subordinator Then Theorem 3.1 
gives the infinitesimal generator of X for t > 0 as 

Afix)=fdaf'ix) + il-f3c)^rix) 

+ /3c / ifix + y)- fix)-yf'ix)lDiy))t^zidy) 

7r\{0} 

+ 7 / [ ifiw) - fix))Pite~‘^,t,xe~‘^'^,dw)izi;idu). 

^ J ( 0 ,cx)) J — 00 

4.4. Martingale with marginals with V symmetric 

Let V be an integrable, symmetric random variable (i.e. V = —V) and {Bt)t>o be a 
Brownian motion independent of V. Following and extending [10], page 283, for any k, 
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let Tt = inf{M > 0: \Bu\ = f^} and Zt = Bt^V. Then {Zt)t>o is a Markov martingale such 
that for each t > 0, Zt = t^V. Moreover, Z is /t-self-similar. Indeed, using the Brownian 

(c) 

motion B^ := cB^-^t, we have 

= inf{u > 0 : | = inf{s > 0 : \Bs\= c-H'^} = c^T^-u^t- 

It follows that = cBc-^c^t _i,„ = cBt and {BT^f)t>ei = {c^BTf)t>o- Hence, 

(Z,,)t>o = = {c^Zt\^^. 

Since Z is a Markov process with transition semigroup 

Po.tfix) = J fit'^v)dF{v), t>0, 

Ps,tf{x) = ^(1 + {s/tY)f{{t/sYx) + i(I - {s/tY)f{-{t/sYx), 0<s<t, 

where F is the cumulative distribution function of H, it has infinitesimal generator 


io/(0)=0, 

Ltfix) = '^xf{x) - yJ{x) + t > 0- 

Using Theorem 3.1, we obtain the infinitesimal generator of the mimic X, for t > 0, 
Affix) = j (^xf{x) + i/3/(-x) - i/3/(x)^ 

+ 7 / [ (/(y) -/(a;))^’(te““,t,a;e"'"“,dy)i/(du). 

^ J (0,oo) J —oo 

Proposition 4.5. The predictable quadratic variation of X is 

Y 1 

(X,X)t = 2K / -Xfds, t>0. 

Jo s 

Proof. Let fix)=x^. Since /(j/)P(te““,t,a;e“'‘“,dy) = and 

Affix) = the result follows immediately from Equation (3.1). □ 


Note that (Z, Z)t = 2k f* ^Zf ds. The predictable quadratic variations of X and Z are 
given by the same functional of the process. 
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4.5. Extension to mimicking Brownian martingales 

Now we discuss how we can (and cannot) alter our martingale condition to mimic some 
Brownian related processes, including the martingales associated with the Hermite poly¬ 
nomials and the exponential martingale of Brownian motion. 

Consider the Hermite polynomials which are defined by 

E vF 

— -hn(x) = expiux — /2) Vu,xSK, 
n\ 

n>0 

equivalently, hn{x) = Let 

Hn{x,t)=t'^^^hn{x/Vt) VxG]R,t>0 

and Hn{x,0) = x". Then where B denotes a Brownian motion, is a local 

martingale for every n since Hn{x,t) is space-time harmonic, that is, -I- \ = 0. 

Take n = 2, the process H 2 {Bt,t) = B^ — t is Markovian and 1-self-similar, thus can 
be mimicked using our mimicking scheme with any subordinator that satisfies '0(1) = 1. 

For n > 3, Hn{Bt,t) is ^-self-similar, but it is not Markovian (see [6]). So we are 
not able to mimic this process by a direct application of the method described above. 
However, a slight modification of our construction proves sufficient to achieve our aim. 

Let Xt be a mimic of the Brownian motion as in Section 4.1 with fc = 0, but without 
the requirement that Xt be a martingale. Then we have the following result. 


Proposition 4.6. For each n, the process Hn{Xt,t) has the same marginal distribu¬ 
tions as Hn{Bt,t) and is a martingale if and only if 'ip{n/2) = n/2, or = 

{sjtYi^. 


Proof. It is obvious that Hn{Xt,t) and Hn{Bt,t) have the same marginals. Here we 
prove the martingale condition. The transition function of X is 


P{s,t,x,dy) = J P{rt,t,{t/sY^‘^r^/'^x,dy)Gs,tidr), 


s<t, 


where P is the transition function of B and Gs,t is the distribution of e Writing 

Pt as the natural filtrations of Hn{Xt,t) and since Hn{Bt,t) is a martingale, we have 


E[HYXt,t)\Ps]=E 


Hn{y,t)P{s,t,Xs,dy)\Fs 


= E 


HYit/sy^V/^X,,rt)Gs,ti<ir)\Ps 


Using that Hn{ax,t) = a"iL„(a;, ^) for any a > 0, we then obtain 


E[HYXt,t)\P,] = (t/s)"/2iL„(X„s)E[e-("/2)C,n(v.)]. 
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□ 

Therefore, in order to mimic the process Hn{Bt,t), we can mimic Bt, with the martin¬ 
gale requirement = 5 changed to tp{^) = §, and then apply the function Hn{x,t) 
to the resulting process. It is of interest to ask whether the above trick extends to other 
space-time harmonic functions. In particular, could this enable us to mimic the geomet¬ 
ric Brownian motion exp(i?i — t/2). Unfortunately, this is not the case. In fact, Hn{x,t), 
n> 1, are the only analytic functions for which this trick works. 

Proposition 4.7. Suppose that H{x,t) = n=o there exists r > I such 

that n < oo, in other words, H{x,t) is analytic on the set (—r,r)^. Sup¬ 

pose further that H is space-time harmonic, so that H{Bt,t) is a martingale. Suppose 
that Xt mimics Bt with the martingale requirement replaced with ip(k/2) = k/2 for a 
positive integer k. Then H{Xt,t) has the same marginals as H{Bt,t) and is a martingale 
if and only if H{x, t) = t^^‘^hk{ xfy/i) where hk{x) is the kth Hermite polynomial. 

Proof. For H{Xt,t) to be a martingale, we must have for any x and s<t, 

H{x,s)= J H{y,t)P{s,t,x,dy) = J H{{t/sy/'^r^^‘^x,rt)Gs,ti<^r), 

that is, 'E[H{xu^^'^rI/'^, suRu)] = H{x,s) for any x, s and m > 1. Letting u = e* and 
Vt =e*“'’*, this is equivalent to ¥\H(xVt^‘^,sVt)\ =H{x,s) for all x and s, or 

OO 00 

^ ^ am.nx^s^. 

m,n =0 m,n =0 

Therefore, for all m, n and t ^ Inr, we must have = am,n- Thus, either 

a^,„ = 0 orE[U,"+™/Vl. 

Recall that E[Ij^] = E[exp(At — A^t)] = exp(—t('0(A) — A)) and there is at most one A 
satisfying Now, choose C,t such that E[V'(^^^] = 1 for a fc € N*. Then, for all 

(m, n) such that m -\- 2n ^ k, am,n = 0. Therefore, 

Lfe/2J 

H{x,s)= ^ ak- 2 n,nx'"~'^'^S^. 

n—0 


Furthermore, 


Lfe/2J 

H{cx,s) = c^ ^ ak-2n,nx‘^~^'^{s/c^y' =c^H{x,s/c^). 

n—0 

By Plucihska [16] and Fitzsimmons [7], H{x,l) is the fcth Hermite polynomial. □ 
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Corollary 4.1. Let Xt be any mimie of Bt in the sense of Section 2 but without the 
martingale requirement. Although the process exp(Xt — t/2) has the same marginal dis¬ 
tributions as exp(Bt — t/2), it is not a martingale unless Q = t, in which case X = B. 


Appendix 

Lemma A.l. Let (Yt)t>o be a Markov process with infinitesimal generator At and ct be 
a deterministic, differentiable, increasing function in t with derivative c't ■ Then the time- 
changed process (Yt)t>o ■= t>o is also a Markov process with infinitesimal generator 
At = c'tAc^. Furthermore, if f is in the domain of A, then f is in the domain of A. 


Proof. Let T be the filtration of Y and F be the filtration of Y. For any bounded 
measurable function g, we have, for s <t, 


For t where the function c is strictly increasing, the infinitesimal generator of Yct is 

E[/(L'c„)|Ycj =x]- f{x) Cu - Ct 


Atf{x) =hm- 

uX-t 




u — t 


■ = AcJ{x)c't. 


If Cu = Ct in a small neighbourhood of t, then cj = 0 and 

E[f{Y,,)\Y,,=x]-f{x) 


Atf{x) = lim- 

uXt 


u — t 


= 0 = AcJix)ct. 


□ 


Lemma A.2. Let (Yt)t>o be a Markov process with infinitesimal generator At and ct be 
a deterministic, differentiable function in t with derivative c't and Ct for any t. Then 
the process {Yt)t>o ■= {ctYt)t>o is also a Markov process and has generator 

c' 

Atf{x) =-Ki/c,At-KcJ{x) + —xf'{x), 

where tTc is an operator defined by 'Kcfix) = f{cx). Furthermore, if iTcf is in the domain 
of A for any c and f is differentiable, then f is in the domain of A. 


Proof. Let F be the filtration of Y and F be the filtration of Y . Let ht be a function 
such that ht{x) = ctx. Since ht is one-to-one, a{hu{Yu),u < s) = cr(Yu,u < s). Therefore, 
for any bounded measurable function g, we have 

E[g(yt)l^«] = ng o ht(Yt)\Fs] = E[g o ht{Yt)\Fs] = E[g o ht{Yt)\Ys] = E[g{Yt)\Ys]. 


The infinitesimal generator of Y is 


Atf{x) = lim( E 

uXt Y 


7I'c„/(K)|>t = —X 
Ct 


-^c„f 


Ct 


+ ]/iu-t) 
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Lemma A.3. Suppose {Zt)t>o is a k- self-similar Markov process. Suppose P{s,t,x,dy) 
and Lt are, respectively, the transition function and infinitesimal generator of Z. Let Zt = 
e~*'^Zgt. Then (Zt)teR is a time-homogeneous Markov process with transition semigroup 

Ptf{x) = J f{y)P{e~*,l,xe~*^,dy) 

and infinitesimal generator 

Lf{x) = Lif{x) - Kxf'{x). 

Furthermore, if f is in the domain of L and differentiable, then it is in the domain of L. 

Proof. By the scaling property of P, we have 

F(Zt e dy\Zs = x)= P(e", e‘, xe"'‘, e‘'‘ dy) = P{eA^-‘^\ l,xeA*-‘^')>^, dy). 

It follows that Z is time-homogeneous and Ptf{x) = //(y)P(e“*, 1,xe“‘'‘, dj/). The gen¬ 
erator of Z can be obtained by applying Lemma A.l and Lemma A.2, and seeing that 
TTeti^e*LgtTTg-tK = Li from Equation (1.1) with c = e* and s = 1. 

Note that for all c > 0, tTc/ is in the domain of L by the scaling property of Z. Thus, 
writing L as the generator of Zt := Z^t, -Kef is in the domain of L by Lemma A.l. □ 

Lemma A.4. Suppose {xt)t>o is a, time-homogeneous Markov process with semigroup 
Pt and generator L, and {r]t)t>o is a subordinator independent of x with drift (3 and Levy 
measure v. Set Yt = Xvt ■ Then the process {Yt)t>e is a time-homogeneous Markov process 
with generator A where 

Af{x)=l3Lf{x)-\- [ {Puf{x)-f{x))n{du). 

J (0,C3o) 

Furthermore, if f is in the domain of L, then it is in the domain of A. 

If rj has zero drift, then Y is a pure jump process. 

Proof. See Sato [17], Theorem 32.1. □ 
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